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We investigate electron transport through single conjugated molecules - including benzenedithiol, oligo- 
phenylene-ethynylenes of different lengths, and a ferrocene- containing molecule sandwiched between two gold 
electrodes with different contact structures - by using a single-particle Green function method combined with 
density functional theory calculation. We focus on the effect of the basis set in the ab initio calculation. It is 
shown that the position of the Fermi energy in the transport gap is sensitive to the molecule-lead charge transfer 
which is affected by the size of basis set. This can dramatically change, by orders of magnitude, the conductance 
for long molecules, though the effect is only minor for short ones. A resonance around the Fermi energy tends 
to pin the position of the Fermi energy and suppress this effect. The result is discussed in comparison with 
experimental data. 

PACS numbers: 73.40.Cg, 72.10.-d, 85.65.+h 



I. INTRODUCTION 

Electron transport through single conjugated or- 
ganic molecules has attracted considerable experimental 
attention ^i^i-^1^1^1^1^1^1^1 ^ ^i because of the convenience of the 
chemical assembly and their delocalized electronic states and 
the resulting small HOMO-LUMO gap which is very useful 
in molecular electronics especially for the operation under 
low bias voltages. Theoretically, to understand the transport 
properties of junctions bridged by these molecules and to 
be able to calculate their conductance are important in view 
of both applications and fundamental physics. Regarding 
the comparison between theory and experiment, however, 
large discrepancies have been found for some conjugated 
molecules between theoretical predictions^^'^"^'^^'^^'^^ and 
experimental conductances ! Many efforts have been 

made on the theoretical sid o^^i^^-^^-^^ ^^^^^^ to understand 
this discrepancy and several possible reasons have been 
discussed. For example, the self-interaction error (SIE) and 
the underestimation of the transport gap in density functional 
theory (DFT) based ab initio calculations,^^'^"^* and the neglect 
of the dynamical correlation effect in the non-interacting 
Landauer formalism.^^'^^ 

In this paper, we investigate a fundamental technical issue 
in all ab initio transport calculations: the effect of the basis set. 
We find that this effect can also cause a significant change, up 
to orders of magnitude, in the resulting conductance, depend- 
ing on the molecular feature and the molecule-lead coupling 
strength. It is shown that for conjugated organic molecules 
with a small HOMO-LUMO gap the position of the Fermi 
energy (Ef) of the leads in the transport gap, which is deter- 
mined by the molecule-lead charge transfer, is sensitive to the 
size of basis set. For long molecules different positions of Ep 
can lead to significant change in the conductance, though for 
short molecules and strong coupling the effect may be only 
minor. We further show, by considering different molecule- 
lead contact structures, that a resonance around the Fermi en- 
ergy tends to pin the position of Ef and thus suppress the 
basis set effect, indicating that this effect will be much less 



important for the strong coupling limit which gives a large 
equilibrium conductance. 



II. COMPUTATION 

The conjugated molecules considered in our calcula- 
tion include a short benzenedithiol (BDT) molecule, oligo- 
phenylene-ethynylenes (OPE) of different lengths, and a 
ferrocene-containing (EC) molecule, sandwiched between 
two gold electrodes, as shown in Figs. 1 and 5. We 
consider two models for the lead, a thin Au(OOl) wire 
as shown in Fig. 1 and a 5x5 Au(lll) surface as 
shown in Fig. 3; Both models are possible to occur in 
experiments^2i2S42&2i28,29,30,3i,32,33,34 ^^^^ consider two dif- 
ferent molecule-lead contacts, through a flat Au surface or 
through an apex Au atom, both of which are also possible to 
occur in experiments . ^i^^i^^i^^i^^i-^-^i-^" ^ 

We adopt a single-particle Green function method^^'^^ 
combining a Landauer formula with ab initio DET^^ elec- 
tronic structure calculation to investigate the electron trans- 
port through the junctions. ^^'^^'^^ In practice, we divide 
an infinitely long Au-molecule-Au system into three parts: 
left lead, right lead, and device region which contains the 
molecule and large parts of the Au leads (see, for instance. 
Fig. 1) for accommodating fully the molecule-electrode in- 
teraction. The self-consistent Kohn-Sham Hamiltonian of the 
device region and the self-energies from the two semi-infinite 
Au leads (the two subsystems are treated on the same foot- 
ing) are used to construct a single-particle Green function 
from which the transmission coefficient at any energy is cal- 
culated. The conductance, G, then follows from a Landauer- 
type relation. The detailed computational techniques have 
been described previously. For the DFT electronic structure 
calculation, we make use of the optimized Troullier-Martins 
pseudopotentials^^ for the atomic cores. The PEE version of 
the generalized gradient approximation (GGA)"^^ is used for 
the electron exchange and correlation. The wavef unctions are 
expanded using a numerical basis set^^ of different sizes: sin- 
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FIG. 1: Structures of two molecular junctions: (a) Au(001)-OPE5-Au(001) and (b) Au(001)-FC-Au(001). Only the extended molecule (device 
region) is shown, which contains 8 atomic layers of the lead surfaces. The lead is a thin a/2 x \f2 Au(OOl) wire. 



gle zeta plus polarization (SZP), double zeta plus polarization 
(DZP), triple zeta plus double polarization (TZDP), and triple- 
zeta plus triple polarization (TZTP). The atomic structure of 
the junctions, including the molecular structure and the first 
two layers of the lead surface, and the molecule- Au separa- 
tion, are fully optimized by minimizing the atomic forces on 
the atoms to be smaller than 0.02 eV/A using the DZP basis 
set. Possible differences from this optimized configuration in 
experiments are not considered. 



III. RESULTS AND DISCUSSION 

To avoid a computational cost which is too large for the 
large size of basis set, we first consider the thin \pi x \pi 
Au(OOl) nanowire model for the lead, as shown in Fig. 1. 
Eight atomic layers of the lead surface are included in 
the extended molecule to ensure good convergence for the 
molecule-lead interaction. The transmission function and the 
equilibrium conductance for the two junctions bridged by the 
OPES and EC molecules are given in Figs. 2 (a) and (b), 
respectively.^^ As can be seen, for both systems varying the 
size of the basis set changes the position Ep the lead 
in the molecular HOMO-LUMO gap. This position is deter- 
mined by the molecule-lead charge transfer; For long-chain 
molecules the charge transfer will be around the contact re- 
gion and form a contact dipole. Different sizes of basis set 
will affect the amount of this molecule-lead charge transfer, 
until it reaches a converged value. For usual ground- state en- 
ergy a small change in the charge transfer will have only a 
minor effect. However, for molecular conductance, this can 
have a significant effect because it depends critically on the 
position oiEp 'm the transport gap, which is very sensitive to 
the charge transferred - the charge in the tail of the broadened 
HOMO or LUMO state. 

For the OPES junction, the SZP basis set leads to enter- 
ing into the LUMO state and, therefore, a very large equilib- 
rium conductance. The DZP basis set begins to put Ep near 
the tail of the LUMO state and reduces the conductance by 
about a factor of 40. When we further go to the TZDP and 
TZTP basis sets, Ep goes further away from the LUMO state 
and the conductance decreases further. The almost identical 
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FIG. 2: Transmission functions of the two molecular junctions 
shown in Fig. 1: (a) Au(001)-OPE5-Au(001) and (b) Au(001)-FC- 
Au(OOl). Different basis sets and the resulting equilibrium conduc- 
tance (in units of GQ=2e^lh) are indicated in the legend. 



results for the TZDP and TZTP basis sets indicate that they 
have reached the converged limit. Qualitatively, the DZP ba- 
sis set already gives a good result though it is not good enough 
quantitatively. 

The result for the EC molecule [Fig. 2 (b)] shows the same 
trend: The SZP basis set leads io Ep entering into the broad- 
ened LUMO state and a large conductance while the DZP and 
TZDP basis sets put Ep m the middle of the gap and a much 
smaller conductance. From SZP to TZDP (TZTP) the relative 
change in conductance is more than two orders of magnitude 
- an effect which has not be found for any other property. 

Next we investigate the basis set effect for conjugated 
molecular bridges with different lengths and also for a more 
realistic Au(lll)-surface electrode. We consider two large 
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FIG. 3 : Structure of the Au( 1 1 1 )-BDT- Au( 111) junction where the 
lead is a 5x5 Au(lll) surface. Only the extended molecule (the 
device region) is shown, which contains 4 atomic layers of each 
lead surface. The structures of the other junctions bridged by 0PE3, 
0PE5, and FC molecules are similar and are not shown here. 



5x5 Au(lll) periodic surfaces as the leads bridged by the 
following molecules of different lengths: BDT, 0PE3, 0PE5, 
and the FC molecule. To reduce the computational cost here 
we consider only the SZP and DZP basis sets and include 
only four atomic layers of the Au(l 1 1) surface in the extended 
molecule, as shown in Fig. 3. To check the convergence with 
the basis set we also carry out a DZDP calculation for the 
0PE3 system. 

From the calculated transmission functions shown in Fig. 4 
one can see the same effect as we discussed previously: the 
two different sizes of basis set affect the molecule-lead charge 
transfer and cause a shift of about 0.5 eV for the position of 
Ef in the HOMO-LUMO gap. The use of the larger DZDP 
basis set in the test calculation for the 0PE3 system leads to 
only a further ~ O.leV shift [see Fig. 4 (b)] and the change in 
the equilibrium conductance is minor, indicating that the DZP 
results are almost converged. For the BDT, 0PE3, and 0PE5 
series, the equilibrium conductance decreases along with the 
increasing molecule length. For the longest one, 0PE5, the 
basis set effect causes about one order of magnitude differ- 
ence in the conductance while for the shorter one, 0PE3, the 
difference is only a factor of 2, and for the shortest one, BDT, 
it becomes minor. This behavior is a result of the observation 
that along with the increase of the molecule length the smaller 
SZP basis set tends to put Ep closer and closer to the LUMO 
state while the larger DZP basis set tends to put it closer and 
closer to the middle of the gap. As a result, the difference 
between the two basis sets grows as the molecule gets longer 

Note that one can see an approximate exponential decay of 
the conductance as the molecule gets longer: from BDT to 
0PE3, and from 0PE3 to 0PE5, the conductance decreases 
approximately by a factor of 20. This is a natural behav- 
ior of coherent tunneling and the decay ratio, /3, is a prop- 
erty of the molecule. Here the (3 value is about 0.21 
which is very close to an available experimental datum, 0.22 
A~^, for a conjugated molecule, carotenoid polyenes.^ The 
agreement between theory and experiment for the exponen- 
tial decay in the case studied here is similar to that found 
in other cases. For instance, experimental^^ and theoretical'^^ 
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FIG. 4: Transmission functions of the four junctions with the leads 
consisting of 5x5 Au(lll) surfaces: (a) Au(lll)-BDT-Au(lll), 
(b) Au(lll)-0PE3-Au(lll), (c) Au(lll)-0PE5-Au(lll), and (d) 
Au(lll)-FC-Au(lll). Basis sets used and the resulting equilibrium 
conductance (in units of Go=2e^/h) are indicated in the legend. 



values of P for another conjugated molecule, oligophenylene, 
are found to be 0.35-0.5 A~^. Also, for alkane and pep- 
tide chains, (3 is in the range 0.7-0.9 experimentally^^ 
and theoretically. ^^'^^ Naturally, the conductance of conju- 
gated oligomers decays more slowly with chain length than 
that of saturated molecules; thus, the values of (3 for conju- 
gated chains are considerably smaller than those for saturated 
chains. 

In the case of the FC molecule the situation is different, 
where the SZP basis set puts the Ef slightly closer to the 
LUMO while the DZP puts it slightly closer to the HOMO, 
and as a result of this symmetry, the two resulting conduc- 
tances become quite close. Thus in general, the strong system 
dependence of the basis set effect indicates that it can cause 
unpredictable numerical errors which are specific to the sys- 
tem investigated. However, this effect can be expected to be 
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FIG. 5: Structures of the junctions bridged by (a)-(d) BDT, (e) 
OPES, (f) OPES, and (g) EC molecules, where the lead is a 5x5 
Au(lll) surface and the molecule-lead connection is through an 
apex Au atom. Only the molecule and the first atomic layer of 
the lead surface are shown. Eor the BDT molecular bridge, differ- 
ent metastable configurations are shown in (b)-(d) with the energy 
difference from the straight connection, (a), listed. Eor the OPES, 
0PE5, and EC molecular bridges, only the most stable configuration 
is shown, which is similar to the BDT case, shown in (d). 



small for short molecules and strong coupling, like for the Au- 
BDT-Au system, where the strong coupling allows lead states 
tunneling into the molecule to meet in the middle (sometimes 
called metal-induced gap states), and as a result, the transmis- 
sion coefficient is fairly large and flat in the gap, as is shown 
in Fig. 4(a). 

In the following we consider a possible experimental sit- 
uation in which the molecule-gold connection is through an 
apex gold atom, as suggested by a STM experiment.^'^^'^^ We 
consider two possibilities: (1) the connection is straight with 
an optimized Au-S bond length, as shown in Fig. 5 (a) for a 
BDT molecule, which is possible to occur when the STM tip 
is pulled from the gold surface, as suggested in Refs.^ HI, 
and |26|; (2) the connection and the molecule-lead separation 
are optimized, as shown in Figs. 5 (d)-(g). The structure opti- 
mization calculation for Au(l 1 1)-Au-BDT-Au-Au(l 1 1) junc- 
tion shows three meta-stable structures, as shown in Figs. 5 
(b), (c), and (d), which lower the total energy with respect to 
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EIG. 6: Transmission functions of the four junctions with the 
leads consisting of 5x5 Au(lll) surfaces, where the molecule-lead 
connection is through an apex Au atom: (a) Au(lll)-au-BDT-au- 
Au(lll), (b) Au(lll)-au-0PE3-au-Au(lll), (c) Au(lll)-au-0PE5- 
au-Au(lll), and (d) Au(lll)-au-EC-au-Au(lll). Results for both 
the straight connection and the optimized connection, as given in 
Fig. 5, are shown. Basis sets used are indicated in the legend. Note 
the resonance peak around the Fermi energy for the BDT, OPES, and 
0PE5 molecule series. 



the straight connection in Fig. 5(a). The most stable struc- 
ture is given in Fig. 5(d). We further find that for the longer 
molecules, OPES, 0PE5, and EC, the most stable optimized 
contact structure is similar to that in Fig. 5(d). Therefore, we 
only give these optimized structures while omiting those with 
a straight connection in Fig. 5 for the longer molecules. 

The transmission functions for the four junctions are given 
in Fig. 6. There are several features to notice in the re- 
sult. First, the presence of the apex gold atom at the contacts 
changes totally the transport property of the junctions. The 
very narrow peaks indicates that the molecule-lead coupling 
is significantly weakened as compared to the case without the 
apex gold atom. This will significantly reduce the current 
for large bias voltages. Second, the relaxation of the con- 
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FIG. 7: Transmission function of the Au(001)-Au-BDT-Au-Au(001) 
junction (Its structure is shown in the inset) given by a Hartree-Fock 
calculation, as discussed in the text. 

tact structure has only a small effect once the connection is 
through an apex gold atom. Third, a sharp resonance peak 
occurs around the Fermi energy except for the FC molecule 
where the peak near Ep is very small [see Fig. 6 (d)]. 

We call this resonance peak the new HOMO state here- 
after. This resonance peak was found previously also by other 
calculations for straight-connected BDT systems. ^^'^^^^^ Al- 
though the present calculation shows that it is a quite gen- 
eral behavior for conjugated molecules connected by an apex 
gold atom, it is unclear whether or not the occurrence of this 
resonance peak has anything to do with issues in DFT, like 
self-interaction error or the absence of the discontinuity in the 
exchange-correlation (xc) potential. To check if it is an ar- 
tifact due to the issues in DFT we perform a test transport 
calculation for a Au(001)-Au-BDT-Au-Au(001) junction (see 
the inset of Fig. 7) using Hartree-Fock (HF) theory which is 
self-interaction free and does not have the issue of xc discon- 
tinuity. The calculation is carried out using a 6-3 IG* basis set 
for H, C, and S atoms and a CRENBS basis set for Au atom, 
and the result is given in Fig. 7. As can be seen, the resonance 
peak around the Fermi energy remains, indicating that it is a 
general behavior instead of an artifact. 

Finally, the basis set effect is significantly reduced com- 
pared to the previous cases without the apex gold atom. The 
small effect of the basis set can be understood by considering 
the entering of the Fermi energy into the new HOMO state. 
Because of it, a small change in the molecule-lead charge 
transfer cannot shift the position of the Fermi energy a large 
amount due to the total charge conserving. This behavior indi- 
cates that the effect of the basis set will be much less important 
in the strong coupling limit which gives a large equilibrium 
conductance. 



IV. COMPARISON TO EXPERIMENTAL RESULTS 

Finally, we would like to make some comparison between 
the present results and experimental data. The present result 
of equilibrium conductance for the Au(lll)-BDT-Au(lll) 
junction (~ 0.1 Go) is too large compared to experimental 
reports. Similar results were also obtained by other DFT-GF 



calculations. ^^'^^'^^ This overestimation is still an open prob- 
lem and can be due to several possibilities: the difference in 
the atomic structure of the junction between the real exper- 
imental situation and the calculation; the self-interaction er- 
ror; and the underestimation of the transport gap in the DFT 
electronic structure calculation 32i Because of the latter issue 
the conductance of the OPES molecule (~ O.OO6G0) and the 
OPES molecule (~ 0.0003 Go) are also expected to be overes- 
timated. 

On the other hand, in real experiments, when a molecule is 
long enough incoherent or inelastic effects may play a role or 
even take over,^^ and, as a result, the conductance will decay 
more slowly or even linearly with increasing molecule length. 
In such case the theoretical calculation for coherent tunnel- 
ing may underestimate the conductance. In a certain range of 
molecule length these two errors may cancel with each other. 

For OPES molecule the present result of conductance is still 
too large compared to some experimental reports ^^'^^ while 
the very small conductance for OPES is in agreement with 
the experimental observation in Ref.O. However, our re- 
sult of a very small conductance (~ 0.00002Go) for the FC 
molecule differs qualitatively from the experimental result^ 
which shows a good conduction.^ ^ So far it is not clear how 
to understand this discrepancy. One possible reason is the 
structure difference between the experiment and calculation. 
For example, in the experiment the molecule-lead connection 
might be through an apex gold atom which will significantly 
increase the current for small bias voltages, as shown by our 
calculation [see Fig. 6 (d)]. Additionally, the scissor mode 
made by the ferrocene moiety (i.e., rotation of the two five- 
member rings with respect to each other) might also allow 
a higher conductance state to be realized, although this might 
not account for the big qualitative discrepancy. ^ ^ Another pos- 
sibility is the chemical absorption on the lead surface in the 
experiments^ which was partially done in solution first. The 
chemical absorption may change the work function of the sur- 
face of the electrode or form a dipole layer on the surface, 
both of which can change the lineup between the Fermi en- 
ergy and the HOMO or LUMO state and therefore cause a 
large increase in the conductance. 



V. SUMMARY 

We have investigated electron transport through differ- 
ent single conjugated molecules of different length by using 
single-particle Green function method combined with density 
functional theory calculation. It is shown that the convergence 
with respect to the size of basis set is important for ab initio 
transport calculations since the position of the Fermi energy 
in the transport gap is sensitive to the molecule-lead charge 
transfer which is affected by different sizes of basis set. It is 
shown that this basis set effect can be dramatic, up to orders 
of magnitude, while it is only minor for short molecules and 
strong coupling. 

Our calculation also shows that a resonance around the 
Fermi energy tends to pin the position of the Fermi energy and 
suppress this effect, indicating that it will be much less impor- 
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tant for the strong coupling limit which gives a large equilib- 
rium conductance. Comparing to experiments, our result for 
the OPE molecules is in agreement with some experimental 
reports. The result for the FC molecule is not in qualitative 



agreement and the reason needs further experimental and the- 
oretical work to clarify. 

This work was supported in part by the NSF (DMR- 
0506953). 
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